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A-optimal  block  designs  for  comparing  test  treatments  with  a  control 

JOHN  STUFKEN* 

Department  of  Mathematics,  Statistics,  and  Computer  Science 
University  of  Illinois  at  Chicago,  Chicago,  Illinois  60680 

Abstract 

We  consider  the  problem  of  comparing  test  treatments  with  a  control  in  a  proper 
block  design.  We  derive  sufficient  conditions  for  the  A-optiinality  of  both  R-type  and 
S-type  designs,  and  demonstrate  how  these  conditions  can  be  used  to  obtain  families  of 
optimal  designs.  We  give  an  example  for  the  construction  of  the  desired  S-type  designs.  A 
table  with  optimal  R-type  designs  (3  <  k  <  10,  k  <  t;  <  30)  is  also  given. 

l.  Introduction 

In  this  paper  we  study  the  problem  of  comparing  a  set  of  test  treatments  with  a 
control  under  the  assumption  that  the  experimental  units  can  be  arranged  in  a  proper 
block  design.  For  a  brief  review  on  the  available  literature  in  this  area  we  refer  the  reader 
to  the  introduction  of  Hedayat  and  Majumdar  (1985).  We  will  use  the  notation  D0[v,  b,k) 
for  the  collection  of  all  connected  block  designs  with  &  blocks  of  size  k  each,  and  based  on  v 
test  treatments  (labeled  1,. . .  ,t>)  and  a  control  (labeled  0).  An  observation  Yijt,  obtained 
by  applying  treatment  t  (0  <  i  <  v)  on  the  experimental  unit  in  block  j  (l  <  j  <  b) 
in  plot  t  (1  <  l  <  A:),  will  be  assumed  to  follow  the  usual  additive  linear  model  without 
interactions: 

Yijt  =  /*  +  *»  +  Pj  +  Ujl' 

We  refer  to  n  as  the  general  mean,  to  r,-  as  the  effect  by  treatment  i  and  to  (3j  as  the  effect 
by  block  j.  The  error  terms  tiji  are  assumed  to  be  uncorrelated,  and  have  a  common 
mean  0  and  a  common  variance  a2.  The  objective  of  the  experiment  is  to  estimate  the  test 
treatment- control  contrasts  T{  —  Jo,  1  <  i  <  v.  We  assume  that  this  is  done  by  using  their 
least  square  estimates,  which  we  will  denote  by  f,-  —  f0.  These  estimates  will  obviously 
be  design  dependent  and  the  problem  that  remains  is  thus,  for  givQn  r,  6  and  k ,  to  select 
a  design  d  (E  Do(v,b,k)  that,  in  some  sense,  gives  us  good  estimates  for  the  contrasts 
of  interest.  In  this  paper  we  will  judge  the  performance  of  a  design  by  the  so-called  A- 
criterion.  A  design  d*  6  Do{v,b,k)  is  called  A-optimal  (in  its  class)  if  it  minimizes 

Va Td(?i  ~  fo) 

%=  l 

over  all  designs  d  €  Do(v,b,k).  Majumdar  and  Notz  (1983)  obtained  a  sufficient  condition 
for  a  design  to  be  A-optimal.  Utilizing  this  condition  we  will  search  for  families  of  A- 
optimal  designs,  both  R-type  designs  (Section  2)  and  S-type  designs  (Section  3). 


*  Research  is  sponsored  by  Grant  AFOSR  85-0320. 


2.  A-OPTIMAL  R-TYPE  DESIGNS 

The  following  definition  first  appeared  in  Bechhofer  and  Tamhane  (1981). 

DEFINITION  2.1.  A  design  d  G  Do(v,b,k)  is  called  a  Balanced  Test  treatment  Incomplete 
Block  Design  (BTIBD)  if  the  following  conditions  are  satisfied: 

1.  d  is  incomplete. 

2.  There  are  constants  Ao  and  Aj  such  that  5Z>=i  ndojndij  =  Ao,  1  <  *  <  v,  and 

—  1  ndi\jndii j  ~  Ai,  1^*1^  *2  ^ 

The  number  0  <  i  <  v,  1  <  j  <  b,  is  the  number  of  times  that  treatment  i  appears 
in  block  j  under  design  d. 

Such  a  design  possesses  a  great  amount  of  symmetry  with  respect  to  the  test  treatments, 
and  it  can  indeed  be  shown  that  Var^f*  —  fo)  is  independent  of  t,  1  <  i  <  v.  Moreover  the 
information  matrix  for  estimating  the  test  treatment  control  contrasts  under  such  a  design 
is  a  completely  symmetric  matrix.  (See  Bechhofer  and  Tamhane  (1981)).  This  symmetry 
was  utilized  by  Majumdar  and  Notz  (1983)  to  obtain  their  main  result.  Before  stating  it 
we  need  one  more  definition. 

DEFINITION  2.2.  By  a  BTIB(v,b,k;t,s)  we  denote  a  design  d  G  D0{v,b,k )  with  the 
following  properties: 

1.  d  is  a  BTIBD. 

2.  There  are  s  blocks  in  d  with  t  +  1  replications  of  the  control,  while  the  remaining  b-s 
blocks  contain  t  replications  of  the  control. 

3.  d  is  binary  in  the  test  treatments. 

THEOREM  2.1.  (Majumdar-Notz  (1983)).  If  3  <  k  <  v  then  a  BTIB(v,b,k;t,s)  is  A- 
optimal  if 

g(t,s)  =  min{g(x,z)  :  (x,z)  G  A}, 

where 

g(x,z)  :=  (v  -  l)2[bvk(k  -  1)  -  (6x  +  *)(v&  -  v  +  k)  +  bx2  +  2 xz  +  z))-1 

+  ((6x  +  z)k  —  (bx2  +  2  xz  +  z))-1, 

and 

A  :=  {(x,z)  :  x  €  {0, 1, . . . ,  [fc/2]  —  l},z  €  {0,1,...,  6}  and  z  >  0  if  x  =  0}. 

(]•]  denotes  the  largest  integer  function). 

In  the  remainder  of  this  paper  we  will  assume  that  the  condition  3  <  k  <  v  is 
satisfied.  We  will  also  adopt  the  terminology  from  Hedayat  and  Majumdar  (1984),  where 
a  BTIB(v,b,k;t,s)  is  called  a  Rectangular-type  (R-type)  design  if  s  =  0  or  6  and  a  Step- 
type  (S-type)  design  if  s  G  {1,. . .  ,6  —  l}.  Through  a  more  detailed  study  of  the  function 
g  in  Theorem  2.1  Hedayat  and  Majumdar  (1985)  obtained  the  following  result  on  optimal 
R-type  designs. 


Theorem  2.2.  A  BTIB(v,b,k;  1,0)  is  A-optimal  if 

(k  -  2)2  +  1  <  v  <  (k  -  l)2. 

The  main  result  of  this  section  is  the  following  generalization  of  Theorem  2.2. 

THEOREM  2.3.  A  BTIB(vfb,k]t,0)  is  A-optimal  if 

(2.1)  (k  -  t  -  l)2  +  1  <  t2v  <  (fc  -  t)2. 

PROOF:  By  Theorem  2.1  it  suffices  to  show  that  (2.1)  implies  that 

(2.2)  g(t,  0)  =  min (g(x,z)  :  (x,z)  6  A). 

From  Hedayat  and  Majumdar  (1985)  it  follows  that  (2.2)  holds  if  the  following  two  condi¬ 
tions  are  satisfied. 

(2.3)  b(k  -  t){{k  -2 1-  l)(v(ifc  -  1)  -  t)2  -  at2(p  -  2t  -  1)) 

<  v(k  -2 1-  l)(p  -2 1-  l)(k  -  1  +  (v  -  2)i) 

and 

(2.4)  b(k  -  t) (at2(p  -  2t  -h  1)  -  (k  -  2t  +  l)(v(k  -  1)  -  t)2) 

<  v{k  -2 1  +  l)(p  -  2t  +  l)(ik  -  1  +  (v  -  2)t), 

where  a  =  (u  —  l)2  and  p  —  v(k  —  l)  +  k.  Since  the  right  hand  sides  in  (2.3)  and  (2.4)  are 
non-negative  these  inequalities  are  obviously  satisfied  if 

(2.5)  (k  —  2t  —  l)(v(fc  -  1)  -  02  ~  «*2(P  -  2t  -  1)  <  0 
and 

(2.6)  at2{p  -2t  +  l)-{k-2t  +  l)(v(Jk  -  l)  -  02  <  0. 

The  left  hand  side  in  (2.5)  equals  —vq2(v),  while  that  of  (2.6)  equals  vqi(v)y  where 

?i(v)  :=  ( k  -  l)tV  -  ((fc  —  t  -  l)(fc  -‘t-t-  l)(fc  -  2t  -I-  1)  +  2 t2(k  -  l))u 
+  2t(lfc  -  t  -  l)(Jfc  -  2t  +  1)  +  t2{k  -  1) 

and 

?2(v)  ’•=  (k  -  l)t2v2  -  ((k  -  t  -  l)(fc  +  t  -  l)(fc  -  2f  -  1)  +  2f2(fc  -  l))v 
+  2t(k  -  t  -  l)(k  -2t-l)+  t2(k  -  1). 

Hence  (2.5)  and  (2.6)  are  equivalent  to  qj(v)  >  0  and  gi(v)  <  0.  Since  both  q\  and  q2  are 
convex  functions  in  v,  to  establish  the  validity  of  Theorem  2.3  it  suffices  to  show  that 
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(0  «i(i)<0 
(ii)  9a(l)  <  0 
(in)  «,((* -im<« 

(iv)  «,(((*  -  <  -  l)2  +  l)/<2)  >  0. 
A  simple  evaluation  gives  us  that 


<7i(  1)  =  —(k  —  t  —  1  )2(*  -  2t  +  1)  <  0, 


which  shows  (i).  Abo 


<72(1)  =  -{k-t-  l)2(fc  -  2t  -  1)  <  0. 

Thb  shows  (ii).  Further  we  obtain  that 

qi((k  -  t)2/t2)  =  ((—2 1  +  1)1 k3  +  (6t2  -  l)*2 

-  ( t 4  +  6 13  +  3 12  -  2t)k  +  2 15  +  3£4  -  t2)/t2. 

Since  k  >  2t  and  since  the  right  hand  side  in  the  above  equality,  viewed  as  a  function  of 
fc,  is  decreasing  on  (2t,oo)  it  follows  by  replacing  k  by  2t  that 


1)3<0, 


establishing  (iii). 
Finally 


»(((*  -  t  -  l)2  +  1  )/J2)  =  ((*  -  21  -  l)((i  -  l)2  -  1*  + l2)  +  *  -  1  )/!2. 

To  see  that  this  is  non-negative  we  merely  notice  that  k  <  v  <  (Jfc  —  t)2/t2  implies  that 
k  >t2  +  t.  Thb  shows  (iv)  and  completes  the  proof  of  Theorem  2.3. 

□ 

We  would  like  to  point  out  that  it  can  be  shown  that  for  fixed  k  and  t  at  most  one 
value  of  v  which  satisfies  (2.3)  and  (2.4)  for  all  b  >  v  is  not  covered  by  the  relation  (2.1). 

The  result  of  Theorem  2.3  enables  us  also  to  obtain  large  families  of  A-optimal  R- 
type  designs  from  families  of  BIB  designs.  We  formalize  thb  in  Corollary  2.1  and  give  an 
example  of  its  use  in  Corollary  2.2.  Other  examples  can  be  found  in  Stufken  (1986a). 

COROLLARY  2.1.  If  there  are  a  BIBD(v,b,k)  and  an  integer  t  £  {1,2, ...,  [k/2]  -  1}  such 
that 

(i)  v  >  k  +  t,  and 

(ii)  (k  -  l)3  +  1  <  t2v  <  k2, 

then  an  A-optimal  BTIB(v,b,k  +  t;t,0)  exists. 

The  proof  is  a  simple  verification  of  (2.1).  For  parameters  satisfying  the  conditions  in 
this  corollary  the  result  gives  a  formal  justification  for  an  idea  of  Cox  (1958),  who  first 
recommended  the  use  of  R-type  designs. 


COROLLARY  2.2.  If  A  is  a  perfect  square  and  v  >  k+ A1/2,  then  asymmetric  BIBD(v,k,  X) 
can  be  embedded  in  a  A-optimal  BTlB(v,v,k  ■+■  A1/2;  A1//2,0J. 


PROOF:  We  verify  the  conditions  in  Corollary  2.1.  It  is  easy  to  verify  that  A1/2  <  [it/2] 
under  the  assumptions  in  Corollary  2.2.  Condition  (i)  is  satisfied  by  assumption,  while  (ii) 
can  be  seen  as  follows: 

Av  =  A(v  -  1)  +  A  =  k{k  -  1)  +  A  <  it2, 

and 

Av  >  k(k  —  1)  >  (k  —  l)2  +  1. 

□ 

We  conclude  this  section  by  giving  a  table  with  all  A-optimal  BTIB(v,b,k;t,0)  designs, 
3  <  k  <  10,  k  <  v  <  30,  whose  A-optimality  follows  from  Theorem  2.3.  The  given  values 
of  6  in  the  table  may  be  replaced  by  any  multiple  of  the  given  number.  The  BIB(v,b,k-t) 
needed  for  the  construction  of  the  given  BTIB  exists  for  all  parameters.  Most  of  them,  or 
their  complements,  can  be  found  in  Table  1.1  in  Hall  (1986). 


3.  A-OPTIMAL  S-type  designs 


r-V- 


Families  of  A-optimal  S-type  designs  are,  with  one  exception  in  Cheng  et  al.  (1986), 
not  available  in  the  literature.  The  reason  for  this  is  twofold.  Partly  it  is  due  to  the  fact 
that  a  characterization  of  families  of  parameters  corresponding  to  S-type  designs  whose 
A-optimality  can  be  concluded  from  Theorem  2.1  is  harder  than  the  analogue  for  optimal 
R-type  designs.  In  addition,  once  a  family  of  such  parameters  has  been  established,  we 
will  still  have  to  answer  the  often  difficult  combinatorial  question  regarding  their  existence 
by  giving  a  method  for  their  construction.  In  the  case  of  R-type  designs  we  could  refer 
to  the  available  literature  on  BIB  designs,  for  S-type  designs  such  an  easy  way  out  is  not 
available.  Cheng  et  al.  (1986)  proved  a  useful  result  for  the  determination  of  families  of 
parameters  of  A-optimal  S-type  designs.  A  formulation  of  their  result  is: 

THEOREM  3.1.  With  g(x ,  z)  and  A  as  in  Theorem  2.1,  if  for  some  t  G  (0,1,...,  [/c/2]  - 1}, 

{1,2 . 6-1} 


1 

(3.1) 

■> 

then 

rSr 

g(t,s)  <  min{g(t,  s  —  l),p(M  +  1)} 


g(t,s)  =  min{g{x,z)  :  (x,z)  €  A}. 


--  /• 
L 


This  result  can  be  used  if  we  determine  parameters  for  which  (3.1)  is  satisfied  and  for 
which  we  can  then  show  the  existence  of  the  corresponding  S-type  design.  Cheng  et  al. 
(1986)  showed  proceeding  in  this  manner  that  an  A-optimal  BTIB(A:2  —  l,q(fc  +  2)(fc2  - 
l),  A;0,7(/c  +  l)(fc2  -  1))  exists  if  A  is  a  prime  or  power  of  a  prime,  where  is  any  positive 
integer.  In  this  section  we  will  derive  families  of  parameters  which  satisfy  (3.1)  for  the 
case  1  =  0.  They  will  include  the  above  result.  We  will  then  conclude  this  section  with  a 
patchwork  construction  for  designs  with  the  obtained  parameters.  Let  us  start  by  looking 
at  5(0,2),  z  €  (0,6].  It  is  not  hard  to  show  that  this  function  attains  its  minimum  either 
at  6  or,  assuming  v  >  4,  at 


(3.2)  zq  :=  bk((v  -  l)(v  +  1)1//2  -  (v  +  l))/(v  +  l)(v  -  3) 

if  Zq  <  b. 

We  are  interested  in  the  latter  case.  It  can  be  verified  that  zo  <  b  if  and  only  if 
v  >  (k  -  l)2  +  1.  Thus  if  this  condition  holds  we  define  s  =  [20]  or  s  =  [20]  +  1  depending 
on  which  of  the  two  gives  a  smaller  value  for  5(0,2).  With  this  choice  of  s,  condition  (3.1) 
will  be  satisfied,  implying  that  a  BTIB(v,b,k;0,s)  is  A-optimal,  if  it  exists.  This  observation 
is  useful  if  a  particular  choice  of  v,  b  and  k  is  under  consideration.  It  is  then  easy  to  obtain 
s  and  an  attempt  can  be  made  to  construct  the  desired  optimal  S-type  design.  A  useful 
first  step  in  such  an  attempt  is  to  verify  whether  the  necessary  conditions  for  the  existence 
of  a  BTIB(v,b,k;0,s)  are  satisfied.  These  are  (see  also  Hedayat  and  Majumdar  (1981)): 


(3.3) 


(k  -  l)s  =  O(mod  t;) 


(3.4)  k(b  -  a)  =  O(mod  v) 

(3.5)  ( k  -  2)(k  -  1  )s/v  4  (k  —  1  )k(b  —  s)/v  =  O(mod  t;  —  1). 

If  one  or  more  of  these  conditions  are  violated,  or  the  non-existence  of  the  desired 
design  is  obtained  through  some  other  argument,  considerations  as  in  Cheng  et  al.  (1986) 
may  lead  to  an  A-optimal  or  highly  efficient  design  in  Do(v,b,k). 

For  the  purpose  of  deriving  families  of  optimal  S-type  designs  the  form  of  s  is  rather 
complicated.  The  problem  is  made  more  tractable  if  we  make  the  additional  assumption 
that  zo  is  an  integer,  and  thus  s  —  zq.  A  design  that  will  be  obtained  under  this  assumption 
has  the  extra  feature  that  any  number  of  copies  of  it  forms  again  on  A-optimal  design.  It  is 
unknown  whether  this  is  true  in  general.  The  following  result  characterizes  all  parameters 
of  interest  for  which  z0  is  an  integer  in  {1, ...  ,b  —  1}. 

T HEOREM  3.2.  The  parameters  for  which  zq  6  {l, . . . ,  b  —  1}  and  for  which  the  necessary 
conditions  (3.3),  (3.4)  and  (3.5)  are  satisfied  are  given  by 

v  =  a2  -  1,  k  =  0,b  =  7«(a  4-  2 )(a2  —  1  )/6e, 

where  a  >  (3  >  3,  8  =  gcd(/3,  a(a  +  2)),  c  =  gcd(06~l ,a2  —  l)  and  7  is  an  arbitrary 
positive  integer. 

PROOF:  For  zq  to  be  an  integer  v  4-  1  must  be  a  perfect  square,  say  v  =  a2  —  1. 
Let  k  =  0  >  3.  For  zo  <  b  it  must  be  that  v  >  (k  -  l)2  4  1,  or  a  >  0.  From 

(3.2)  we  obtain  now  z0  =  b0(a  4-  l)/a(a  4  2).  Since  gcd(a  4  l,o:(a  +  2))  =  1  we  see 

that  this  is  only  an  integer  if  b  =  O(mod  a(a  +  2)/ 6).  If  we  set  b  =  -f1a(a  4  2 )/6 
we  obtain  with  s  =  zo  from  (3.3)  and  (3.4)  (which  imply  kb  -  s  =  O(mod  v))  that 
7!/?(a2  4  a  —  1  )/6  =  0(mod  a2  —  1),  or  equivalently  ii(3/6  =  0(mod  a2  -  l).  There¬ 
fore  71  =  7(a2  -  l)/e,  for  some  positive  integer  7.  This  gives  the  parametrization  for 
b  as  asserted  in  the  theorem.  It  can  indeed  easily  be  verified  that  (3.3),  (3.4)  and  (3.5) 
are  now  satisfied,  while  zq  =  s  =  7/?(or  4  l)(c*2  —  which  is  indeed  an  integer. 

□ 

This  result  tells  us  thus  that  a  BTIB(a2  -  l,7a(a42)(a2  —  l)/5e,/?;0, 7/?(a4  l)(a2  - 
1  )/6e)  is  A-optimal  for  any  a,  0,  7,  6  and  e  as  in  Theorem  3.2.  Although  the  conditions 

(3.3) ,  (3.4)  and  (3.5)  are  also  satisfied,  this  does  unfortunately  not  guarantee  the  existence 
of  the  design.  Patchwork  constructions  may  resolve  this  problem,  either  for  individual 
designs  or  for  certain  families  of  designs  of  the  above  form.  We  give  here  one  example 
of  such  a  construction.  Let  a  =  0  =  k.  This  results  in  a  BTIB(&2  -  l,7(fc  4  2)(&2  - 
1) ,  A;  0, 7 (k  4  1  )(k2  -  1)),  the  existence  of  which  was  shown  for  k  a  prime  or  power  of  a 
prime,  as  referred  to  earlier.  We  add  to  this  the  following  result. 

THEOREM  3.3.  An  A-optimal  BTIB(k 2  -  l,~j(k  +  2)(k2  -  l),&;0,7(/:4  l)(£2  -  1))  exists 
if  k  4  1  is  a  prime  or  power  of  a  prime. 


PROOF:  It  suffices  to  show  the  existence  for  7  =  1.  The  layout  of  the  desired  design  is  as 
follows: 


1 


i 


1 

K 


(k+1) (k  -1) 


1|  control 


We  have  to  determine  the  parts  labeled  by  A  and  B.  Start  by  partitioning  the  k2  —  1  test 
treatments  into  k- 1  groups  Gi, ...,  Gk-i  of  cardinality  A  +  l  each.  For  each  of  these  groups 
form  the  trivial  BIBD,  (fc-t-l)Efc.  This  gives  a  total  of  k2  —  1  blocks  of  size  k,  blocks  that  we 
use  to  form  B.  Next  form  an  orthogonal  array  OA{[k+\)2,  k— l,fc+l,  2)  of  index  unity.  The 
construction  of  this  array  is  well  known  if  k  + 1  is  a  prime  or  power  of  a  prime,  which  is  just 
our  assumption.  Use  this  orthogonal  array  to  construct  a  group  divisible  design  with  k2  -  l 
treatments  in  blocks  of  size  k  —  1  in  which  the  treatments  are  partitioned  in  k  —  1  groups  and 
a  pair  of  treatments  from  the  same  group  does  not  appear  at  all,  while  treatments  from  dis¬ 
tinct  groups  appear  exactly  once  as  a  pair.  Clearly  we  can  arrange  this  such  that  the  groups 
are  again  G 1, . . . ,  G*-i.  Take  k  —  1  replications  of  this  group  divisible  design  and  let  this  be 
A.  It  is  easy  to  verify  that  the  design  constructed  in  this  way  is  the  desired  S-type  design. 

□ 

Several  other  patchwork  constructions  for  individual  designs  can  be  found  in  Stufken 
(1986a). 

4.  Concluding  remarks 

The  results  in  Section  2  extend  the  main  results  of  Hedayat  and  Majumdar  (1985). 
Theorem  2.3  and  Corollary  2.1  provide  a  sufficient  condition  for  the  A-optimality  of  R-type 
designs,  a  condition  that  can  easily  be  verified.  Large  families  of  optimal  R-type  designs 
can  be  obtained  from  it.  It  is  unfortunately  known  that  for  many  families  of  parameters 
v,  b  and  k  the  best  design  is  not  a  R-type  design.  An  alternative  would  then  be  to  search 
among  S-type  designs.  Section  3  gives  some  families  of  optimal  S-type  designs  and  some 
ideas  for  finding  optimal  designs  for  specified  parameters.  The  construction  of  these  S- 
type  designs  is  usually  a  non-trivial  problem.  Patching  BIB  designs  and/or  group  divisible 
designs  together  in  an  appropriate  way  may  resolve  problems  of  this  kind.  After  these 
considerations  there  will  still  be  many  parameters  v,  6,  and  k  for  which  no  optimal  design 
has  been  found.  This  will  be  the  case  if  optimal  designs  are  neither  of  the  R-type  nor 
S-type  or  if  Theorem  2.1  is  not  strong  enough  to  determine  their  optimality.  In  such  cases 
one  could  either  search  for  highly  efficient  R-type  designs  (see  Stufken  (1986b))  or  give 
consideration  to  the  ideas  suggested  in  Cheng  et  al.  (1986). 


We  finally  point  out  that  all  the  optimal  designs  derived  in  this  paper  are  not  only 
A-optimal,  but  also  MV-optimal. 
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